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Abstract 

The cosmological significance of the neutralino sector is studied for a class of super- 
symmetric grand unified theories in which electroweak symmetry breaking is seeded 
by a gauge singlet. Extensive use is made of the renormalization group equations 
to significantly reduce the parameter space, by deriving analytic expressions for all 
the supersymmetry-breaking couplings in terms of the universal gaugino mass 
the universal scalar mass tuq and the coupling A. The composition of the lightest 
supersymmetric partner is determined exactly below the W mass, no approxima- 
tions are made for sfermion masses, and all particle exchanges are considered in 
calculating the annihilation cross-section; the relic abundance is then obtained by 
an analytic approximation. We find that in these models, stable neutralinos may 
make a significant contribution to the dark matter in the universe. 



(to appear in Nuclear Physics B) 



1 Introduction 



The neutralino (x) in supersymmetric theories is a leading candidate for the constituent 
of the dark matter in the universe [|]. It has been demonstrated to be the lightest 
supersymmetric partner (LSP), and to have a significant cosmological relic density, in both 
the minimal supersymmetric standard model (MSSM) P|-|[L0||, and in an extended version 



TT[| in which symmetry breaking is driven by a gauge singlet [|12|-|14[]. In contrast to other 
dark matter candidates such as the axion or gravitino, the neutralino has a relic density p x 
such that its contribution to the cosmological density parameter, Q x = p x j p c , is naturally 
close to unity [|TJ. Here p c ~ 1.05 x lCT 5 /i 2 GeV cm -3 is the critical density corresponding 
to a flat universe and h = H /100 km s _1 Mpc -1 is the (conventionally scaled) Hubble 
parameter, which is observationally constrained to lie in the range 0A<h<, 1.0 

The actual value of Q is rather uncertain. Luminous matter, e.g. in the disks of 
galaxies, contributes at most Q ~ 0.01 while studies of dark matter in galactic halos and 
in groups and clusters of galaxies imply Q ~ 0.01 — 0.2 [|TIJ. Indirect estimates of Q based 
on observations of peculiar (i.e. non-Hubble) motions on supercluster and larger scales 
suggest higher values in the range ~ 0.2 — 1 [|HJ. The recent observations by the COBE 
satellite of angular fluctuations in the microwave background radiation has given strong 
support to the hypothesis that the universe is dominated by dark matter with a density 
close to the critical value; interpretations of these data require that for cold dark matter, 
Q ~ (0.2 — 0.5) h~ l flI8| . Additionally, the product Qh 2 (oc p) can be conservatively 



bounded from above by requiring that the age of the universe exceed the observational 
lower limit of ~ 10 10 yr; this yields Vtti 1 <, 1 for h > 0.4, assuming a Friedmann-Robertson- 



Walker cosmology [15]. Further, if one assumes that the dark matter provides the critical 
density, then the age constraint restricts the Hubble parameter to be h<2/3, hence 
flh 2 ,$,4/9. Taking all this into account, we adopt as the range of cosmological interest: 

0.1 < a x h 2 <0.5 . (1) 

Note that the range 0.01 <, fl x h 2 <, 0.1 would be relevant only to the dark matter in 
galactic halos. Obviously the constituent of the cosmological (cold) dark matter may also 
be the constituent of the halo dark matter (but not neccessarily vice versa). 

In the standard picture of the evolution of the universe |15j, it is assumed that ex- 



pansion is adiabiatic following an early period of inflation. Massive particles 'freeze out' 
of chemical equilibrium with the thermal radiation-dominated plasma when their rate of 
annihilation falls behind the rate of expansion; hence their relic abundance is inversely 
proportional to the thermally-averaged annihilation cross-section (cr ann t>). In order to ex- 
pedite the calculation of this quantity in supersymmetric theories, most authors resort 
to assumptions about the values of various parameters in order to reduce the multi- 
dimensional parameter space (e.g. 20-dim in the case of the MSSM !). One common 
assumption, motivated by grand unified theories (GUT), is that the gaugino masses are 
related by their renormalisation from a common value at an energy scale of O(10 16 ) GeV, 
viz. 
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m Y = - Mr )m 2 . (2) 



One can go further and note that all the low energy couplings are then related by the 
renormalization group equations (RGEs) and a few model-dependent parameters. Use 
of this fact is however seldom made in the literature, because the RGEs are not soluble 
analytically, and also because the large value suggested by experiment for the mass of 
the top quark requires its Yukawa coupling to be quite close to the unitarity limit in the 
MSSM, thus prohibiting any acceptable approximations. 

On the other hand, attention has recently been drawn H [| |TU| to the importance of 



using accurate expressions for the low energy couplings when determining the relic abun- 
dance. The model considered in these papers was the MSSM as derived from supergravity- 
inspired GUTs, which has only 6 arbitrary parameters defined at the unification scale. 
These are the universal scalar (m ) and gaugino (mi/ 2 ) masses, the Higgs coupling [fi) 
and the hidden sector parameters (A and B) and the top quark mass (m t ). In prin- 
ciple, the remaining parameters such as the Higgs vacuum expectation values (VEVs) 
are determined by the electroweak breaking constraints, although in practice matters are 
complicated by large radiative corrections to the effective potential from top and stop 
diagrams [Tj|. 

In this paper we shall carry out a similar analysis for the 'minimally' extended super- 
symmetric standard model. Running the RGEs analytically, we have obtained polynomial 
approximations for all the parameters at the weak scale in terms of their values at the 
GUT scale. Comparison with numerical evaluations show our solutions to be accurate to 
better than 5%, even close to the unitarity limit. This reduces the dimensionality of the 
parameter space to 6: m^j, m Q , A, the top quark Yukawa coupling A2 t , and two Yukawa 
couplings involving the Higgs fields and gauge singlets, A7 and A 8 . Again the Higgs VEVs 
are determinable from the electroweak symmetry breaking requirements. In the following 
section we describe the model and introduce the dark matter candidate. Then we review 
the prescription for analytically obtaining the relic density to an accuracy better than 5% 
and, using our RGE approximations, calculate its value for various choices of the Higgs 
VEVs. We find that the neutralino has a relic abundance Q x ~ 0(1) over most of the 
experimentally allowed parameter space. 

2 The Minimally Extended Standard Model 

The model we shall consider is the MSSM extended by including a gauge singlet to 



give the necessary Higgs mixing and thus break SU(2)l ® U(1)y 0, 01 • The relevant 
superpotential is 

W = A 1 D R (D[h°-U[h-) + A 2 U R (U[h°-D[h + ) (3) 

+A 3 E R (E[h° - i/Jh") + A 7 (h°h° - h-h + )0 o + A 8 4 (4) 
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The first three terms give masses to the (s)quarks and (s)leptons, such that Ai(h), Aa(h) 
and As(h) are the mass matrices of the down quarks, up quarks, and electrons respectively 
(the generation indices have been supressed). The form of eq. (|3]) is clearly similar to 
the Standard Model superpotential, apart from the presence of the last two couplings 
involving the gauge singlet 0o, which replace the usual quadratic Higgs coupling, and 
ensure a satisfactory breaking of SU(2)l <S> U(1)y- This model is the minimal extension 
to the MSSM that contains only trilinear couplings, thus making it a natural candidate 
to emerge from level-one fermionic string theory. (It should be pointed out that in the 
'flipped string' models |22] no gauge singlet appears in the low-energy effective potential, 
because such singlets may aquire large VEVs at a scale of C(10 n ) GeV corresponding to 
the condensation of scalars in the hidden SU(4) ® SO(10) sector. We are therefore not 
addressing schemes of this form.) 

The phenomenology of such models has been studied extensively in the literature 
|12|]-|[14]], [plH-|p7ll, but relatively little attention has been paid to the cosmological relic 
abundance of the neutralinos therein. Other proposed dark matter candidates are the 
'flatino' superpartner of the massless scalar flat field ('flaton') after SU(5) <8> U(l) sym- 
metry breaking [EBI and the integer-charged bound states ('cryptons') in some hidden 



sector [p7 ]. The flatino/neutralino mass matrix is restricted by the amount of R-parity 
breaking which is induced. This is because models with broken R-parity (such as the min- 
imal flipped SU(5) model considered in ref. fl26f ) have relic particles which are unstable. 
The decays of such particles into photons, charged particles or even neutrinos may violate 
the astrophysical bounds obtained by analysis of the diffuse 7-ray background and data 
from underground and cosmic ray detectors |28| . Careful analysis of the diagonalisation 



of the neutral fermion mass matrix shows |29[] that the mixing of the neutralino with the 
neutrino is 0(mw/m G u T ) 2 , not the generic value of 0(m w /m G u T ) adopted earlier |26fl . 
Assuming that flaton decay repopulates the neutralino density, we find that the neutralino 
lifetime is ~ 10 6 — 10 15 yr and falls in the region forbidden by the astrophysical bounds 
We are therefore forced to reject the flatino as a dark matter candidate, at least in 



the form presented in ref. ||26|| . There are few phenomenological restrictions on cryptons at 
present since their nominal cosmological abundance would be far in excess of cosmological 
limits and must needs be substantially diluted by invoking entropy generation subsequent 
to freeze-out. Hence their relic density is quite uncertain | |2^| . Further, such particles 
can also decay through non-renormalizable superpotential interactions and are thus con- 
strained by the aforementioned astrophysical bounds if their relic density happens to be 
comparable to that of the dark matter |28[ . 

Inspired by the most general supergravity models []30[ , we shall break supersymmetry 
softly at a scale mx with the initial conditions of universal gaugino masses mi/2, scalar 
masses mo and trilinear scalar couplings A. For the sake of expediency we specifically 
assume m x = hjgut- However, we argue that changing the particular details of the 
model at high energies has only minor consequences for our analysis, and that we are 
therefore addressing a large class of possible theories. Consider, for example, the case 
where mx = lCT 2 mGUT- Between mx and mcuT, the gauginos and squarks receive 
contributions from the RGEs of ~ 5% of their final values, so that the assumed degeneracy 
of the soft breaking parameters is approximately correct. Likewise, we expect the effect 
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of the deconfinement of any hidden sector to be mitigated if it happens reasonably close 
to the GUT scale. In addition, we observe the following: any parameters which do not 
depend on the values of the strong coupling or top quark Yukawa coupling, converge 
to infra-red stable points at low energy and are therefore relatively independent of the 
(dimensionless) details of the model at high energies. On the other hand parameters 
such as A 7 which do depend on the strong coupling or top quark Yukawa coupling, are 
dominated by their (diverging) values at low energies and are thus also expected to be 
relatively stable to changes in the details of the model at high energy. 

We shall only consider interactions relevant for neutralino masses below that of the 
W boson. (Above the W mass, the annihilation cross-section receives important contri- 
butions from processes such as XX ~~ > W + W~ [|J). Since the infra-red fixed point of As is 
only ~ 0.21 [11], we expect our analysis to be valid for most of the parameter space where 
(0o) an d mi/2 are not too large. When the neutralino mass exceeds mw, our results are 
only qualitative. 

If we neglect the bottom quark mass, then the squark and slepton mass matrices which 
we need are already diagonal, and follows: 

the up and charm squarks, 



/ ™2 



m 



m 2 Q -g 2 (v 2 -v 2 )/12 + g 2 (v 2 -v 2 )/4 

m 2 v + g\{y 2 - v 2 )/3 
the down, strange and bottom squarks, 



UL 



\ ur 



(5) 



/ m 2 - \ 



the left and right handed sleptons, 



g 2 (v 2 -v 2 )/12-g 2 (v 2 -v 2 )/4 
m D — 9y(v 2 — v 2 )/6 



777 



cl 



V ml 



and, the sneutrino, 



ml + g Y (v 2 - v 2 )/A - g\{v 2 - v 2 )/A 



mi 



-g 2 Y (v 2 -v 2 )/2 



(6) 



(7) 



"4 = mi + g Y (v 2 - v 2 )/A + g 2 (v 2 - v 2 )/A 



(8) 



where v, v and x are the VEVs of the scalar fields h°, h° and 6n, and the rest of the 



notation is as in ref. . 

In addition we have three Higgs states, whose mass squared matrix is: 
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( g 2 v 2 -(A 7 X 7 +3X 7 X 8 x)vx/v (2X%-g 2 )vv+(A 7 +3X 8 x)X 7 x (A 7 +6X 8 x)X 7 v+2X%vx\ 

(2X 2 --g 2 )vv+(A 7 +3X 8 x)X 7 x g 2 v 2 - (A 7 X 7 +3X 7 X 8 x)vx/v (A 7 +6X 8 x)X 7 v+2Xjvx 

\(A 7 +6X 8 x)X 7 v+2X 2 vx (A 7 +6X 8 x)X 7 v+2X 7 vx 3A 8 X 8 x+36Xjx 2 -A 7 X 7 vv/x J 

(9) 

where g 2 = (g 2 , + g 2 )/2, and we define, for future use: tan/3 x = x/v, tan/5 = v/v. Finally 
there are two physical states coming from the pseudoscalar mass squared matrix, 



/ (A 7 + 3X 8 x)X 7 xV 2 /vv (A 7 X 7 - 6X 7 X 8 x)V 

\ (A 7 X 7 - 6X 7 X 8 x)V (A 7 + 12X 8 x)X 7 vv/x + 9A 8 X 8 x 



(10) 



where V 2 = (v 2 + v 2 ). The A 7 and A 8 are the coefficients of the A7 and X 8 trilinear scalar 
coupling terms. In order to be able to identify the Goldstone boson state in the above, we 



have substituted in the minimisation conditions |T]J , and find the massless eigenvector to 
be 

G = sin ph. - cos (3h. , (11) 

as expected from consideration of current conservation. This use of the minimisation 
condition has the added advantage that it enables one to avoid estimating the SUSY 
breaking, gauge singlet mass-squared term, which is the most sensitive to small changes 
in the Yukawa couplings. As stated earlier, the remaining low energy parameters may be 
approximated to better than 5% using the RGEs as follows. 

We begin at m^y and renormalize the Yukawa and gauge couplings to mguT where c/3 
and #2 are unified. We have used 2-loop RGEs for the gauge couplings, and 1-loop RGEs 
for all other parameters, as in ref. p3| . At each step in the renormalization procedure, 
we keep track of the parameters by expanding them in terms of the Yukawa couplings at 
rn-yv- Then the supersymmetry breaking terms mi/2, A and m are switched on, and all 
the parameters evolved back down to m w . This gives a set of polynomials in Ai, A 2 , A 7 , 
As, ^1/2, mo and A which we truncate at sixth order in the Yukawa couplings. 

The current LEP values of the strong coupling (extracted from Z event shapes using 
resummed QCD) |32| and of the weak mixing angle |31| are given by : 

a s (m z ) = 0.124 ±0.005 , sin 2 6 w (m z ) = 0.2325 ± 0.0007 , (12) 

where both quantities are evaluated in the MS scheme, as is appropriate for our renor- 
malization group calculations. However, analyses of deep inelastic /i, v scattering give a 
somewhat smaller value: a s (mz) = 0.112±0.005, while analyses of r decays and of J/\I/, T 



decays give intermediate values, all with comparable errors [32]. Given the present uncer- 



tainty in determinations of a s , we prefer to perform our calculations using the 'weighted 



mean' value 32 
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a s {m z ) = 0.118, (13) 

with sin 2 9w(mz) = 0.2325, and present the corresponding results in Appendix A. (Note 
that the equivalent expressions for the MSSM may be deduced by setting A7 = As = 0, 
replacing A7 by B^, and the first term in by Brno, where B is the quadratic scalar 
Higgs coupling at the GUT scale.) In order to illustrate the effect of a small change in 
a s , we also solve the RGEs using a s — 1/8 = 0.125 in Appendix B and using a s = 0.113 
in Appendix C. It should be noted that since v/|<7y is still significantly below g 2 at the 
GUT scale, the value of my is slight grey larger than suggested by eq. (£3). In fact the 
correct relation is: 



-»v= S (14) 




These approximations turn out to be much better than expected, again because of the 
convergence of the couplings to infra-red stable points. Thus, despite the fact that the 
analytic approximations diverge from the correct values during the running of the RGEs 
towards the GUT scale, upon returning to the weak scale we find an accuracy of better 
than 5% when compared with numerical runnings, even close to the unitarity limit. The 
results obtained are sensitive to the value of a s and the current uncertainty (±0.005) in 
its determination corresponds to an uncertainty of ~ 10% in the calculated couplings. 
We find a discrepancy of ~ 10% with the results obtained using one-loop expressions for 
the gauge couplings (which for the first two generations reproduce the numerical values 
obtained in ref. ||). For the case a s = 1/8, we find good agreement with the numerical 
results of ref. 133], taking A 2t = 0.31. 



Relic Abundance of the Neutralino 



The relic neutralino energy density may be computed by solution of the Boltzmann trans- 
port equation governing the abundance of any stable particle in the expanding universe. 
The exercise is greatly simplified by noting that under suitable assumptions this reduces 
to the continuity equation |34] 



dV x 



n, 



71 



cq\2 



R 



3 . 



(15) 



this simply states that the number of neutralinos (of density n x ) in a comoving volume 
is governed by the competition between annihilations and creations, the latter being 
dependent, by detailed balance, on the equilibrium neutralino density, n^ q . This equation 
holds if the annihilating particles are non-relativistic and if both the annihilating particles 
and the annihilation products are maintained in kinetic equilibrium with the background 
thermal plasma through rapid scattering processes |35]. . The essential physical input is 
the 'thermally-averaged' annihilation cross-section, which may be conveniently calculated 
in the rest frame of one of the annihilating particles (i.e. the 'lab' frame) 
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/■oo 

(ffamw) = / de JC(x,e) a ann v , (16) 

"2x s — Arv? 

/C(x,e) = -^yv^(l + 2e) ir 1 (2x v / TTi) , x ee m x /T , e = — x , 

where t> is the relative velocity of the annihilating neutralinos, s is the usual Mandelstam 
variable, and the K n are modified Bessel functions of order n. This 1-dim integral can 
be easily evaluated numerically and is well-behaved even if <7 ann t> varies rapidly with e, 
as near a resonance or at the opening up of a new annihilation channel. (Useful analytic 



approximations have been provided for these situations in refs. [36, p7|.) Outside these 



problematic regions, a simple series expansion for (cr ann f) is useful [B6l. BH]: 



(a,n«v) = a<°> + \a^x^ + ( V + ^a< 2 )) x~* (17) 

+ ['^ a (D + l^ a (2) + ^ a (3)^-3 + . 



16 16 16 

where oW is the nth derivative of (a ann v) with respect to e, evaluated at e = 0. We will 
be interested in the value of (cr ann v) for very non-relativistic annihilating particles (with 
x ~ 20), hence it is adequate to retain terms only upto 0{x~ 1 ) if errors of a few percent 
are tolerable. Therefore we write 

(^^)-a (0) + ^ (1) ^. (18) 

Since the neutralino (x) is a Majorana fermion, (the S-wave part in the above) is 
suppressed relative to aP^ (which contains both S and P-wave pieces) due to the smallness 
of the final state fermion masses 0. As noted above, this expansion breaks down close to 
resonances where the value of a^ 1 ' may even be negative; however careful analysis p6|, |37 



shows that the true situation is not in fact pathological. We are mainly interested in 
identifying the broad characteristics of parameter space, eq. ( |T8"D is therefore adequate for 
our purposes. The above expansion also breaks down close to particle thresholds because 
one of the denominators in decreases to zero. In our case this happens when the x 
mass is close to that of the bottom quark; however as we will see, such low masses for 
m x are in fact disallowed by experiment. Finally, we stress again that only annihilation 
channels relevant for m x < m w are considered here. Above the W mass, (<7 ann v) receives 
important contributions from processes such as XX W + W~, as discussed in ref. ||. 
This is relevant to the present discussion only for small values of tan f3 x and large values 
of mi/2, where our results are only qualitative. 

Most solutions to the continuity equation fllSP in the literature (e.g. |15|, 0, |35| ) 
have been obtained by transforming variables from t —>■ T assuming that the product of 
the cosmic scale-factor R and the photon temperature T is constant. However, RT does 
change as the temperature in the cooling universe falls below various mass thresholds and 
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the corresponding particles annihilate and release entropy (i.e. thermalized interacting 
particles) into the system, according to 



dR 



dT 



3 g ai 



(19) 



Here g Sl counts the number of interacting degrees of freedom which determine the specific 
entropy density, sj: 



g SI (T) 

where, si 



45 Sl (T) 



2tt 2 T 3 



(20) 



3m 2 + 4p 2 
3Ei(p)T 



f? q (P,T) 



d 3 p 



Here the sum is over all interacting particles, Ei = y m 2 + p 2 , gi is the number of inter- 
nal (spin) degrees of freedom, and the equilibrium phase-space distribution function (for 
zero chemical potential) is, f^ q (p,T) = [exp(Ei/T) ± with +/— corresponding to 
Fermi/Bose statistics. 

In addition we have to keep track of any particle species which 'decouples', i.e. goes 
out of kinetic equilibrium with the thermal plasma, while still relativistic, thus taking its 
entropy content out of the 'interacting' sector (which defines the photon temperature). 
The quantity which remains truly invariant (in a comoving volume) is the total entropy 
S = sR 3 , where 



s(T) 



2tT 
45 



9s(T)T 3 



with g s , the entropy degrees of freedom, related to g Sl as p6| 

g s AT Dj 



g s {T) = g Sl (T)U 

dec 



1 + 



(21) 



:22) 



The product is over all the decoupled particle species j (which decouple at T = T Dj ). 
Such particles will not share in any subsequent entropy release hence their temperature 
will drop below that of the interacting particles. By conservation of the total entropy, we 
obtain \\ 
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T, 



T 



g Sj (Tv) g Sl (T) 



1/3 



9 S] (T) g Sl (T D )\ ' ^ 

The total energy density includes, of course, all particles, decoupled or interacting, ap- 
propriately weighted by their respective temperatures. The energy degrees of freedom 
9 P , is given by 



g P {T) 



30p(T) 



7T 



2 J" 4 



(24) 
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where, p = 9k E k (q) fl q (p,T k 

all J 



d 3 p 



In the present case the only particles which decouple while relativistic are massless neu- 
trinos, hence g s is actually the same as g si above the neutrino decoupling temperature of 
O (MeV), and can replace the latter in eq. ([19|) (see e.g. ref. [ j38|] ) 

The continuity equation (P~5| ) can now be rewritten [^] in terms of the dimensionless 
variables Y x = n x /s, Y£ q = nf^/s, and x = m x /T as 



dx 



Ax- 2 [(Y?) 2 - F x 2 j , (25) 
where, A = m x m P (cr ann v) g l J 2 , 

and, ,y 2 (T) ^ J0lf 1 + ldlng SI (T)\ 



9 P {T) 



3 dlnT 



The parameter keeps track of the changing values of both g p (which measures the 
total energy density and thus determines the expansion rate R/R), as well as g Sl (which 
determines the specific entropy density and consequently the value of the adiabat RT). [] 
The values of g p (T) and g Sl (T) have been computed in ref. |38| by explicit integration over 



the particle distribution functions (using eqs. 120 and El), and g*(T) has been calculated 



from these data in ref. [36]. There is considerable ambiguity concerning the behaviour 



of these quantities during the quark-hadron phase transition and the authors of ref. [B£ 



present curves for two choices of the critical temperature, T c , viz. 150 MeV and 400 MeV 
(corresponding to two choices of the 'bag' constant characterizing the confined hadrons). 
In fact the two curves merge at T ■> 800 MeV, above which temperature it is appropriate 
to treat the thermal plasma as an ideal gas of non-interacting fermions and gauge bosons. 
As we shall see, the cosmologically interesting regions of parameter space correspond 
to neutralino masses which are high enough (m x >, 20 GeV) that freeze-out occurs well 
above the quark-hadron phase transition; hence their relic abundance is not particularly 
sensitive to the choice of T c . To simplify the presentation of our results, we specifically 
adopt the curves corresponding to T c = 150 MeV. Another source of ambiguity concerns 
the electroweak phase transition, since if the relevant critical temperature is small then 
this may alter the relic abundance of massive particles which obtain their mass at this 
phase transition |39|] . The LEP bounds on the Higgs mass imply that this does not happen 
in the minimal Standard Model fl39| , however in supersymmetric models, the electroweak 
phase transition may well be strongly supercooled. Nonetheless this ought not to affect the 
relic abundance of neutralinos lighter than a few hundred GeV, which is of interest in this 
work. (However, such effects ought to be taken into account for the heavier neutralinos, 
m x ~ O (TeV), considered in refs. [g § 0.) 

According to the continuity equation ([25]), the particle abundance tracks its equilib- 
rium value as long as the self-annihilation rate is sufficiently rapid relative to the expansion 

1 Warning: is often used (e.g. |l5[| ) to denote the energy degrees of freedom, which is called g p here. 
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rate. When the particle becomes non-relativistic, its equilibrium abundance falls expo- 
nentially due to a Boltzmann factor, hence so does the annihilation rate. Eventually, the 
annihilation rate becomes sufficiently small that the particle abundance can no longer 
track its equilibrium value but becomes constant (in a comoving volume). Hence the 
parameter 

A = ye^ , (26) 

X. 

grows exponentially from zero as the particle goes out of chemical equilibrium, according 

Ahg K 2 {x) K 2 (x) l dln^ 

4^^T A A(A + 2) ~ Kjx) ~ x dhrT ' (27) 

where we have used Y x q = (A5g / An 4 )[x 2 K 2 (x) / g S {T)]. This equation can be solved for the 
'freeze-out' temperature, x = Xf T corresponding to a specific choice of A = A fr ; it is found 



that the choice Af r = 1.5 gives a good match to numerical solutions [5B|. In fact, the 
first term on the r.h.s. is ~ 1 for non-relativistic particles, while the second term on the 
r.h.s. is zero and gl = g~ 1 ^ 2 if g Sl is constant at the epoch of freeze-out; the freeze-out 



temperature is then given by the more commonly used condition (e.g. JT5|) : 



x b = ln[A fr (2 + A fr ) 5] - ^\nx iT (28) 
/ 45 \ 1/2 

where 5 = {—^ ) g m x m P (a ann v) g p 1/2 . 



7f 



This corresponds to the epoch when the annihilation rate, nf q (cr ann v ), equals the rate 



of change of the particle abundance itself f34[: dlnn eq /dt = R/R. We find that for 
neutralinos, Xf r ~ 20 — 25. 

After freeze-out, only annihilations are important since the temperature is now too low 
for the inverse creations to proceed; hence the asymptotic abundance, Y Xoo =Y x (t — » oo), 
obtains by integrating dY x /dx = —\x~ 2 Y x with the initial condition, Y x (xf r ) = Y x q (xf r ) 
33. This gives 



(29) 



45 [ I \ 1/2 

— m P / {cr ann v) g J dT 

7T JO 



Since Y Xoo subsequently remains invariant, the number density of the neutralinos today 
(at T = T ) is obtained by simply multiplying by the present entropy density, s(T ), 
which is given by eq.(|l]) with g s (T ) = 3.91 and T = 2.735 ± 0.02 K g|. The present 
neutralino energy density is thus given by 
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(7TI \ 
GevJ ■ (30) 

The observational error in T corresponds to an uncertainty of ±2% in Q x h 2 . The above 
procedure give an accuracy of better than 5% when compared with numerical solutions 
of eq.(|2|) (e.g. |§). 

Two further remarks are in order. We have assumed above that the neutralino, as the 
LSP, is absolutely stable by virtue of R-parity invariance. However, most models have 
been considered as effective theories derived from broken GUTs; in some cases (e.g the 
minimal flipped SU(5) model) a small amount of R-parity breaking is possible, causing 
the LSP to decay via superheavy boson exchange. As noted earlier, if such metastable 
LSPs constitute the dark matter, then constraints from astrophysical data |28j require 
their lifetime to be at least 10 6 times longer than the age of the universe, hence they may 
be considered effectively stable. 

We have also assumed that the neutralino is much lighter than the charginos and 
sleptons, so that it cannot annihilate into these particles at the epoch of freeze-out; such 
'co-annihilation' effects can greatly reduce the relic abundance ||37|| . In fact such effects are 
important and must be taken into account for pure Higgsino states ref. fl4lj| . However in 
the present case the neutralino is either a gaugino or a mixed state over most of parameter 
space. Also, the presence of the gauge singlet ensures that there is no degeneracy between 
the neutralino and the squark masses as argued in ref. [|]J for the case of the MSSM. 



4 Exploring Parameter Space 

We now take the most general expression for the cross-section and composition of the 
LSP and analyse the parameter space in the (A 7 , A 8 ), (m , mi/ 2 ), (tan/5, tan/5 a .) and 
(tan/?, A) planes. Below the W mass, the calculation of the cross-section is simplified 
since the contributing diagrams are only those shown in fig. 1 (where f denotes a matter 
fermion). The cross section may therefore be found as a general expression involving 
left and right external fields and arbitrary internal propagators; summing over all mass 
eigenstates and interactions yields the required quantity (which is far too cumbersome to 
be presented here). We have checked that our cross sections for various pure LSP states 



agree with those found in refs. [|1J, [3£ 



We restrict our attention to SUSY breaking parameters less than 500 GeV, to avoid 
excessive fine-tuning in the model. As discussed earlier, we look for cosmologically 
'interesting' regions in which the value of the relic density is 0.1 <, Q x h 2 <, 0.5, and exclude 



regions in which Q x h 2 >1. Using recent experimental data from LEP and CDF 



we also exclude regions of parameter space where the slepton, sneutrino, chargino, 
(lightest) Higgs or gluino have masses which have already been excluded (or are less than 
that of the neutralino), or where the supersymmetric contribution || ^5f to the (recently 
revised) Z° 'invisible width' j|5[] is excessive. In particular, we impose the following bounds 
(all at the 90% confidence level): 

2 For a model- independent survey see ref. p2|. 
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m,~ v ± > 45 GeV 



m h > 43 GeV, m 9 > 40 GeV, nig > 150 GeV, 



(N„ -3) < 0.11 . 

(31) 

We also adopt m t = 120 GeV as suggested by recent analyses of high precision LEP data 
46fl . It is neccessary to exclude negative Higgs and pseudo-Higgs masses-squared since 



otherwise the absolute minimum of the potential would be CP-violating, as is implicit in 
their mass matrices f|TT 



4.1 The (A 7 , As) plane 

Since the sfermion masses are independent of A7 and As to first order (see the appendices), 
when we look in the (A 7 , A 8 ) plane we see the effect on the LSP relic density of the masses 
of the LSP, Higgs and pseudo-Higgs only. (In the two remaining planes, however, the 
masses of all the (s) particles in the low energy spectrum have an effect.) For a gaugino- 
like LSP, sfermion exchange is the important contribution to (cr an nv), but Z exchange 
becomes increasingly important for Higgsino-like and mixed LSPs. Also, for Higgsino-like 
LSPs, the sfermion exchange diagrams are suppressed by a factor of at least m b /m z . 

In contrast with previous work, it is not possible here to choose the signs of Aj and 
A 8 , in order to guarantee the positiveness of the Higgs and pseudo-Higgs masses-squared. 
Instead we choose the signs of the VEVs such that the experimental conditions above are 
satisfied, and allow A7 and A 8 to be both positive and negative. This translates to the 
following constraints on the VEVs: 



Sgn[v] = +1 

Sgn[v] = -Sgn[A 7 A 8 ] 

Sgn[x] = +SgnL4A 8 ]. (32) 

We find no other choice of parameters which consistently meet all experimental require- 
ments. The first condition is chosen arbitrarily, without loss of generality. Inspection of 
the neutral terms in the potential shows that under the above conditions, all the masses 
are invariant under the transformations A7 — > — A7 and A 8 — > — A 8 . We are therefore free 
to choose A7 and A 8 to be positive. 

In figs. 2a-d we show contours of the LSP mass and the relic density in the (A 7 — A 8 ) 
plane, as well as the experimentally allowed regions. (Note that A7 and A 8 are bounded 
by their infra-red fixed values of 0.87 and 0.21 respectively). We choose tuq and m\/2 
to be comparable and of order the Fermi scale since the changing LSP mass is then the 
dominant parameter in determining the relic density. For small values of A7 and A 8 , the 
LSP mass is small and the relic abundance large. However these regions are unlikely to 
satisfy experimental bounds and we find that this is indeed so. At large A7 and A 8 , the 
Higgs and pseudo-Higgs masses increase and therefore do not contribute appreciably to 
(cr ann i;), which is then dominated by sfermion exchange. Also the LSP becomes mostly 
gaugino, and its mass becomes large and nearly constant. For positive A, the relic density 
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tends to a plateau (at a relatively low value) in the experimentally allowed region (see 
fig. 2a). (The 'holes' seen in the figures correspond to poles in the annihilation cross- 
section where the calculated relic density is neligibly small. As noted earlier, this is an 
artifact of the approximate series expansion eq. (|T8|); correct thermal averaging of the 
cross-section at resonances shows that the dips in the relic abundance are not in fact as 
sharp or as deep |36|, |37j.) 



The behaviour is similar, although somewhat more complex, for negative A, but the 
more mixed LSP states produce a significantly lower relic density in the experimentally 
allowed region. Higher values of Q x h 2 occur along the strips where A7 ~ 0.6, and where 
A7 ~ 0.25, when the LSP mass is still relatively small (see fig. 2b). Because of the relations 
in eq. (|32"1), the relic density takes longer to reach its plateau at high values of the Yukawa 
couplings since the LSP eigenstate is dominated by the mixed higgsino/singlet. 

In addition, making \A\, tan (3 or m large causes the relic density to flatten out more 
quickly, and tends to increase its value in the flat region (see figs. 2c, d). We find that 
the experimentally allowed region tends to shrink for higher values of tan/3 and \A\ and 
therefore conclude that these are experimentally disfavoured. These last two observations 
will be discussed in more detail below. 

For the values of parameters chosen in fig. 2d (i.e. A = +1, tan/3 = 4, tanf3 x = 1/2, 
m t = 120 GeV, m = 300 GeV and m 1/2 = 200 GeV), fi x h 2 flattens out at 

fi x /i 2 ~0.4, (33) 

as seen in fig. 2e. In the experimentally allowed region, there is very little variation with 
As, for either positive or negative A. In this region the neutralino mass lies in the range 
- 60 -80 GeV. 

Clearly the large value of A7 (= g 2 ) often chosen in the literature |14| is in fact favoured 



by experiment; also the relic density in this region is particularly constant. The value of 
As ~ ^2/6 [[14] however seems to be too low, and Ag ~ 0.2 is a better choice. (In any case 
the former value was derived using smaller values of the top quark coupling, and we find 
that in our case the latter value is quite natural for A2 t of 0(1) at the GUT scale.) 

Henceforth, we shall choose A 7 = g 2 and A 8 = 0.2, to analyse the (m , m^) and (tan /3, 
tan f3 x ) planes, bearing in mind that when the LSP contains a significant proportion of 
higgsino or gauge singlet, the values obtained for Q x h 2 may be slight grey reduced by 
choosing larger A 7 . 



4.2 The (mo, mi/2) plane 

In figs. 3a-d we exhibit contours of the relic density, together with the experimentally 
allowed region, and the mass of the LSP (which depends only on the value of m^). 
There are three easily identifiable regions in the (mo, mi/2) plane. The first is the region 
close to the origin, where the relic density is small. This is because here the sfermion, 
Higgs and pseudo-Higgs masses are small, causing (cr ann t;) to be correspondingly large. 
The LSP mass increases linearly with m^, until it becomes dominated by the gauge 
singlet, when it becomes constant. At this point the relic density tends to increase, 
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as the sfermion exchange contribution to (<J ann v) (the third diagram of fig. 1) becomes 
suppressed. This increase becomes more marked for larger tan/3, and as a result, the 
region of high my 2 ( £ 300 GeV) appears to be excluded for these values. The third 
region one can identify is where mo is large, the LSP mass is small, and the Higgs and 
pseudo-Higgs masses are large. This also yields a cosmologically interesting relic density 
for m £ 200 GeV, which is relatively independent of tan (3 and A. This behaviour is 
similar to that found for the MSSM in ref. ||. At low the LSP is mostly photino 
and has a mass of order | sin 2 6\ym 2 + ©(mj/mz), in accord with ref. JL4J. 

We find that high values of fl x h 2 are favoured by high values of mo and low values of 
mi/2- In particular we note that for tan /3 £ 4 and negative A, there are no allowed regions 
which are cosmologically interesting. For positive A, the no-scale theories (mo = 0) are 
excluded by LEP as shown in fig. 3a. Note that since we are using the choice of VEVs 
given by eq. fl52"D, the sign of A may instead be regarded as a choice of the sign of A 8 . For 
negative A this is not the case, and we find that there exists an upper limit on the allowed 
mo, coming from the requirement that the vacuum not be CP- violating (see fig. 3b). This 
bound becomes more restrictive if the values of tan/5 or \A\ are increased as seen in 
figs. 3c, d. 

4.3 The (tan/3, tan/^) and (tan/3, ^4) planes 

Choosing the above values of A7 and As, we examine the relic density in the (tan /3, tan /3 X ) 
plane (figs. 4a,b) for A = ±1, my 2 = 200 GeV, and m = 300 GeV. The experimental 
constraints discussed earlier now place the following limits on the allowed values of these 
parameters, 

0.5 < tan/3< 5.0 

tan ft < 2.0, [6V 

corresponding to the values of A7 and As chosen here. (Higher values of tan (3 considerably 
shrink the experimentally allowed regions in fig. 2.) These limits are considerably less 



restrictive than those suggested in ref. ||14j| , where the specific values chosen were A7 = 
and As = #2/6 — 0.108. In addition we find that, remarkably enough, the parameter space 
is restricted to those values of tan/3 and tan j3 x which yield f2 x ~ 1, and where the relic 
density is relatively uniform, especially for positive A. In fig. 5, we examine the effect of 
varying A and tan /3. Now Q x h 2 is relatively constant over the whole plane, apart from 
at A = where we find a jump in the relic density from Q x h 2 ~ 0.3 down to Q x h 2 ~ 0.1 
due to the change in the sign of v and v dictated by eq. (|3^). In general we conclude that 
negative values of A favour higher values of fl x . We can place additional limits on the 
trilinear coupling, 

-3<A<4, (35) 
which marginally favour positive values of A. 
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5 Conclusions 



Using a set of approximate analytic solutions to the RGEs, we have examined the relic 
density of the LSP in a class of unified theories in which the electroweak symmetry is 
broken by a gauge singlet. We have used the most general expressions for all the sparticle 
masses and for the thermally-averaged cross-section below the W mass. This is particu- 
larly important when the LSP is gaugino-like, in which case the dominant contribution to 
the cross-section is the sfermion exchange. In addition, we stress the importance of using 
accurate expressions for the scalar and trilinear coupling terms A-j and As, since these 
are involved in determining the electroweak breaking. The unknown Yukawa couplings 
A 7 and A 8 are experimentally restricted to have relatively high values, 

0.2 <A 7 < 0.8 , . 

0.1 <A 8 < 0.21, [6Q) 

where the last limit is a unitarity constraint, above which the naive unification assumption 
is no longer valid. For the chosen values of m t = 120 GeV, A7 = #2, A 8 = 0.2, the values 
of the VEVs are constrained by eq. ([34] ) and the trilinear coupling A at the GUT scale 
by eq. (|3"5]). Further, positive values of A are marginally favoured by experiment and 
lead to a more uniform relic abundance in parameter space. We find that cosmologically 
interesting regions exist with 



250 GeV < m < 500 GeV. (37) 

For negative A, higher values of mo tend to be excluded by symmetry breaking require- 
ments. The region m 1( / 2 £ 300 GeV is excluded for larger tan /3, although one should now 
take into account new processes involving W final states, and would therefore expect fl x h 2 
to be smaller than our estimates in this region. As pointed out in ref. [13], for various 
'pure' LSP states (no gaugino content) there are cosmologically interesting regions for 
much higher values of m^, although such large masses (exceeding 1 TeV) necessitate 
fine tuning of at least one part in a hundred in the top quark Yukawa coupling, to obtain 
satisfactory electroweak breaking ||42|] . (We also confirm the conclusion of ref. [jlj that 



for a higgsino-like LSP (corresponding to large mi/2 and mo and small A 8 ) there is no 
bound on the LSP mass coming from the constraint fl x h 2 ;$ 1.) Although computational 
constraints dictated a tree-level treatment of electroweak symmetry breaking in this work, 



we note that 1-loop corrections have a relatively minor effect on the relic density ||14|| . We 
therefore anticipate that our conclusions will remain unaltered in a more detailed analysis 
of this model including radiative corrections. 
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Appendix A 



Analytic approximations for all required parameters (at mz) in terms of the Yukawa cou- 
plings (at mz) and the supersymmetry-breaking terms m-o, m.1/2 and A for the minimally 
extended, supersymmetric, standard model, with a s = 0.118. Note that ml, Wq, ml, Ai 
and A2 carry suppressed generation indices. Q 

A x = Amo(l-0.78A?-0.13A|-0.31A?+0.03A?A?-0.02AlA?+0.25A?A|) 

+ m 1/2 (-3.85+1.69A?+0.01A^+0.28A^+0.11A?+0.04A?A?-0.02A^A?+0.25A^) 

A 2 = Am (l-0.13A?-0.77Al-0.31A^-0.02A?A^+0.03A^+0.25A?A|) 

+ m 1/2 (-3.89+0.28A?+1.70A^+0.01A^+0.11A?-0.02A?A?+0.04A^A^+0.25A^) 

A 7 = Am (l-0.39A?-0.39A^-1.23A?-0.11A?A2-0.11A^A^-7.46A^-0.52A?A^ 
+ m 1/2 (-0.58+0.85A?+0.85A^+0.43A?-0.10A?A?-0.10A^A?-0.50A^) 

A 8 = Am (l-1.84A?-0.22A^A?-0.22A^A?-22.4A^-3.05A?A| 
+ m 1/2 (0.69Af-0.2lA?A^-0.21A^Af-3.01A^) 

ml = m2(l-0.77A2-0.26A 2 A2+0.03A 2 A2A2+0.20A 2 A^+0.07A2A2+0.13A 2 A^A2) 
+ Am m 1/2 (1.14A^-0.15A 2 A^-0.86A^-0.25A^A^) 
+ m 2 /2 (6.27-3.21A 2 ,+0.15A 2 A 2 .+0.94A 2 1 -0.07A 2 ,A 2 .) 

m 2 Q = m 2 (l-0.39A 2 -0.13A 2 A 2 + 0.10A 2 Af-0.39A 2 -0.13A 2 A 2 + 0.03A 2 A 2 A 2 + 0.10A 2 A^ + 
0.03A 2 A?+0.06A 2 A 2 A 2 +0.04A 2 A 2 +0.06A 2 A 2 A 2 ) 
+ Am m 1/2 (0.56A 2 -0.44Af+0.57A^-0.14A 2 A^-0.43A^-0.13A 2 A?-0.13A^A?) 
+ m 2 /2 (6.68-1.59A 2 +0.46A^-1.60A 2 +0.16A 2 A 2 +0.47A^-0.03A 2 A 2 -0.03A 2 A?) 

ml = m 2 (l-0.78A 2 -0.26A 2 A 2 +0.20A 2 A^+0.03A 2 A 2 A^+0.07A 2 A?+0.13A 2 A 2 A?) 
+ Am m 1/2 (1. 13A 2 - 0.87A? - 0. 14A 2 A 2 -0.25A 2 A 2 ) 
+ m 2 /2 (6.22-3.17A 2 +0.93A^+0.16A 2 A 2 -0.07A 2 A 2 ) 



m 2 E 


= mQ + 0.15m 2 / 2 


ml 


= mQ + 0.51m 2 / 2 


my 


= 0.41m 1/2 


m 2 


= 0.79m 1/2 


m 3 


= 2.77mi /2 



3 The equivalent expressions for the MSSM may be obtained by setting A7 — Ag = 0, replacing A7 by 
and the first term in by Bm , where B is the quadratic scalar Higgs coupling at the GUT scale. 
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Appendix B 



As in Appendix A, with a s = 0.125. 

A 1 = Am (l-0.76A?-0.12A^-0.31A?+0.04A?A?-0.02AlA?+0.26A?A|) 

+ mi /2 (-4.09+1.75A?+0.01Af+0.29A^+0.11A?+0.04A?A?-0.02AlA?+0.26A?Ag) 

A 2 = A?7io(l-0.13A?-0.75A^-0.31A?-0.02A?A?+0.04AiA?+0.26A?Ag) 

+ mi /2 (-4.13+0.29A?+1.75A^+0.01A^+0.11A?-0.02A?A?+0.04AlA?+0.26A?Ag) 

A 7 = A7no(l-0.38A?-0.37A|-1.23A?-0.11A?A?-0.11A^A?-7.51Ag-0.48A^+0.32A|) 

+ mi/ 2 (-0.59+0.87A?+0.88A^+0.44A?-0.11A?A?-0.11AlA?-0.52A?Al) 

A 8 = A77io(l-1.85A?-0.22A?A?-0.22A|A?-22.5A|-3.0A^+0.95A|) 

+ m 1/2 (0.66A^-0.22A?A^-0.22A^A^-3.10A^) 

ml = ml{l-0.75Xl-0.25A 2 Xl+0mA 2 XlXl+0A9A 2 X^+0^XlX 2 7 +0A2A 2 XlX 2 7 ) 

+ Am m 1/2 (1.17A^-0.15A 2 A 2 -0.87A^-0.26A^A^) 

+ m 2 /2 (6.90-3.44A 2 i +0.17A 2 A 2 . + 1.01A 2 1 -0.07A 2 3 A 2 .) 

m 2 Q = m 2 (l-0.38A 2 -0.13A 2 A 2 + 0.10A 2 At-0.37A 2 -0.12A 2 A 2 + 0.03A 2 A 2 A 2 + 0.10A 2 A^ + 
0.04A 2 A?+0.06A 2 A 2 A?+0.04A^A^+0.06A 2 A^A?) 

+ Am m 1/2 (0.58A 2 -0.44Af+0.58A 2 -0.15A 2 A 2 -0.43A^-0.13A 2 A 2 -0.13A 2 A 2 ) 

+ m 2 /2 (7.30-1.71A 2 +0.50A^-1.72A^+0.17A 2 A^+0.50A^-0.03A 2 A 2 -0.03A^A?) 

ml = m 2 (l-0.76A 2 -0.25A 2 A 2 +0.19A 2 A^+0.03A 2 A 2 A^+0.07A 2 A?+0.13A 2 A 2 A?) 

+ Am m 1/2 (1.16A 2 -0.88A?-0.14A 2 A 2 -0.26A 2 A 2 ) 

+ m 2 /2 (6.84-3.41A 2 +1.0Af+0.17A 2 A^-0.07A 2 A^) 



ml 


= ml + 0.15m 


ml 


= ml + 0.52m 


m Y 


= 0.40mi/ 2 


m 2 


= 0.78m!/ 2 


m 3 


= 2.89m 1/2 
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Appendix C 

As in Appendix A, with a s = 0.113. 

A 1 = Am (l-0.80A?-0.13A^-0.30A?+0.03A?A?-0.02AlA?+0.24A?A|) 

+ mi /2 (-3.65+1.64A?+0.01Af+0.28A^+0.10A?+0.04A?A?-0.02AlA?+0.24A?Ag) 

A 2 = A?7io(l-0.13A?-0.79A^-0.30A?-0.02A?A?+0.03AiA?+0.24A?Ag) 

+ mi /2 (-3.69+0.27A?+1.65A^+0.01A^+0.10A?-0.02A?A?+0.04AlA?+0.24A?Ag) 

A 7 = A7no(l-0.40A?-0.39A|-1.21A?-0.11A?A?-0.11A^A?-7.33A|-0.47A^+0.15A|) 

+ mi/ 2 (-0.57+0.82A?+0.83A^+0.42A?-0.10A?A?-0.10AlA?-0.49A?Al) 

A 8 = A77io(l-1.82A?-0.21A?A?-0.21A|A?-22.0A|-2.84A?Al+0.46A|) 

+ m 1/2 (0.62A^-0.21A?A^-0.21A^A^-2.87A^) 

m^j = m2(l-0.79A2-0.26A 2 A2+0.04A 2 A2A2+0.20A 2 A^+0.07A2A2+0.13A 2 A^A2) 

+ Am m 1/2 (1.10A^-0.14A 2 A^-0.86A^-0.24A^A^) 

+ m 2 /2 (5.74-2.99A 2 + 0.15A 2 A^ + 0.89A^-0.06A 2 A^) 

= m 2 (l-0.40A 2 -0.13A 2 A 2 + 0.11A 2 Af-0.39A 2 -0.13A 2 A 2 + 0.03A 2 A 2 A 2 + 0.10A 2 A^ + 

0.03A 2 A?+0.06A 2 A 2 A?+0.03A^A^+0.06A 2 A^A?) 

+ Am m 1/2 (0.55A 2 -0.43Af+0.55A 2 -0.14A 2 A 2 -0.43A^-0.12A 2 A 2 -0.12A 2 A 2 ) 

+ m 2 /2 (6.15-1.48A 2 +0.44A^-1.50A^+0.15A 2 A^+0.44A^-0.03A 2 A 2 -0.03A^A?) 

ml = m 2 (l-0.80A 2 -0.27A 2 A 2 +0.21A 2 A^+0.03A 2 A 2 A^+0.07A 2 A?+0.13A 2 A 2 A?) 

+ Am m 1/2 (1.09A 2 -0.86A?-0.14A 2 A 2 -0.24A 2 A 2 ) 

+ m 2 /2 (5.70-2.97A 2 +0.88A^+0.15A 2 A^-0.06A 2 A 2 ) 





= mQ + 0.15m 2 / 2 


ml 


= mQ + 0.50m 2 / 2 


m Y 


= 0.42mi /2 


m 2 


= 0.79m 1/2 


m 3 


= 2.67m!/ 2 
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Figures 



Figure 1 Diagrams contributing to the thermally-averaged cross- sect ion. 

Figure 2a Contour plots of Q x h 2 and the LSP mass (in GeV) in the A 7 — A 8 plane, 
for A = +1, tan/5 = 2, t&n/3 x = 1/2, m t = 120 GeV, m = 300 GeV and m 1/2 = 
200 GeV. The cosmologically interesting region 0.1 < Q x h 2 < 0.5 is shaded light 
grey while the forbidden region Q x h 2 > 1 is shaded black; the dark grey regions 
correspond to intermediate values of f2 x h 2 . The bold hatched line encloses the 
experimentally allowed region. The arrows on the axes indicate the infra-red fixed 
points. 

Figure 2b As in fig. 2a, but with A = — 1. 
Figure 2c As in fig. 2a, but with A = —3. 
Figure 2d As in fig. 2a, but with tan f3 — 4. 

Figure 2e Q x h 2 versus A 7 with A 8 = 0.2; the other parameters are as in fig. 2d (A = +1, 
tan/3 = 4, tan/3 x = 1/2, m t = 120 GeV, m = 300 GeV and m 1/2 = 200 GeV). The 
experimentally allowed region is demarcated by the vertical hatched lines. 

Figure 3a Contour plots of il x h 2 and the LSP mass (in GeV) in the m — plane, 
for A = +1, tan (5 = 2, tan(3 x = 1/2, A 7 = g 2 , A 8 = 0.2 and m t = 120 GeV. The 
cosmologically interesting region 0.1 < Q x h 2 < 0.5 is shaded light grey while the 
forbidden region Q x h 2 > 1 is shaded black; the dark grey regions correspond to 
intermediate values of Q x h 2 . The bold hatched line encloses the experimentally 
allowed region. 

Figure 3b As in fig. 3a, but with A = — 1. 

Figure 3c As in fig. 3a, but with A = —3. 

Figure 3d As in fig. 3b, but with tan (3 = A. 

Figure 3e Q x h 2 versus mo with mi/2 = 200 GeV; the other parameters are as in fig. 3b 
(A = -1, tan = 2, tan(3 x = 1/2, A 7 = g 2 , A 8 = 0.2 and m t = 120 GeV). The 
experimentally allowed region is demarcated by the vertical hatched lines. 

Figure 4a Contour plots of Q x h 2 and the LSP mass (in GeV) in the tan f3— tan f3 x plane, 
for A = +1, A 7 = g 2 , A 8 = 0.2, m t = 120 GeV, m = 300 GeV and m 1/2 = 200 GeV. 
The cosmologically interesting region 0.1 < VL x h 2 < 0.5 is shaded light grey while 
the forbidden region VL X h 2 > 1 is shaded black; the dark grey regions correspond 
to intermediate values of Vt x h 2 . The bold hatched line encloses the experimentally 
allowed region. 

Figure 4b As in fig. 4a, but with A = — 1. 
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Figure 5 Contour plots of £l x h 2 and the LSP mass (in GeV) in the tan (3 — A plane, 
for tan = 1/2, A 7 = g 2 , A 8 = 0.2, m t = 120 GeV, m = 300 GeV and m 1/2 = 200 
GeV. The cosmologically interesting region 0.1 < Q x h 2 < 0.5 is shaded light grey 
while the forbidden region Q x h 2 > 1 is shaded black; the dark regions correspond 
to intermediate values of Q x h 2 . The bold hatched line encloses the experimentally 
allowed region. 
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